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Abstract 

In these three lectures I discuss the present status of high-energy astroparti- 
cle physics including Ultra-High-Energy Cosmic Rays (UHECR), high-energy 
gamma rays, and neutrinos. The first lecture is devoted to ultra-high-energy 
cosmic rays. After a brief introduction to UHECR I discuss the acceleration of 
charged particles to highest energies in the astrophysical objects, their propa- 
gation in the intergalactic space, recent observational results by the Auger and 
HiRes experiments, anisotropics of UHECR arrival directions, and secondary 
gamma rays produced by UHECR. In the second lecture I review recent results 
on TeV gamma rays. After a short introduction to detection techniques, I dis- 
cuss recent exciting results of the H.E.S.S., MAGIC, and Milagro experiments 
on the point-like and diffuse sources of TeV gamma rays. A special section 
is devoted to the detection of extragalactic magnetic fields with TeV gamma- 
ray measurements. Finally, in the third lecture I discuss Ultra-High-Energy 
(UHE) neutrinos. I review three different UHE neutrino detection techniques 
and show the present status of searches for diffuse neutrino flux and point 
sources of neutrinos. 



1 Ultra-high-energy cosmic rays 
1.1 Introduction 
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Fig. 1: Left: The cosmic ray spectrum I{E) as function of kinetic energy E, compiled using results from the 
LEAP, proton, Akeno, and HiRes experiments 1 1 , 2|. The energy region influenced by the Sun is marked in yellow 
and an 1/ E'^-^ power-law is also shown. Right: The same spectrum at high energies E > 10^^ eV multiplied by 
E^ |[3|. Spectrum changes are called the 'knee' at 10^^ eV and the 'ankle' at 10^^ eV. 



Particles coming from space to the atmosphere of the Earth historically were called cosmic rays. 
Most cosmic rays, however, are not 'rays' or photons, but charged particles, protons and nuclei. Real 
high-energy gamma rays coming from space to the Earth are only a small fraction of total flux, and they 
will be discussed in Section[2] The measured spectrum of cosmic rays from 100 GeV to highest energies 
E > 10^^ eV is presented in Fig. [l] (left). The yellow strip at low energies presents the contribution of 
the Sun. The remaining spectrum can be fitted with a single power law 1/E'^-^ up to highest energies. 



The main contribution to it above 100 GeV gives galactic sources. After multiplication of the spectrum 
on the energy cube, one can see changes of power law in Fig. [l] (right). Ai E > 10^^ eV the spectrum 
becomes steeper. This change in the spectrum called the 'knee' and associated energy E = 10^^ eV 
is the maximum energy up to which galactic sources accelerate cosmic rays. The next change of the 
spectrum is located at £^ = 3 • 10^^ eV and has two possible interpretations. Either this is the place where 
extragalactic sources start to dominate or it is the result of pair- production energy loss by extragalactic 
protons (see Section [OJ. At the end of the spectrum there is a cutoff, which was not seen in the old 
experiments presented in Fig. [T] (right) due to small statistics, but it was observed recently by the HiRes 
Q and Auger ||4j| experiments. 

In this lecture I briefly discuss the theory and observations of Ultra-High Energy Cosmic Rays 
(UHECR), the highest-energy particles measured on Earth with energy E > 10^^ eV. Such particles, 
protons and nuclei, can be accelerated in astrophysical objects, propagate through intergalactic space, 
losing energy in the interactions with Cosmic Microwave Background (CMB). UHECR are charged 
particles. Therefore they are also deflected in the Galactic and intergalactic magnetic fields on the way 
from the source to the Earth. For a more detailed introduction to UHECR I recommend recent lectures 
by M. Kachelriess ||5|. 

There are several important scales commonly used in astroparticle physics. Distance is usually 
measured in parsecs, 1 pc = 3 • 10^^ cm. Corresponding larger units are kiloparsec 1 kpc = 10^ pc and 
megaparsec 1 Mpc =10^ pc. Energy at highest energies is usually expressed in units of EeV = 10^^ eV. 

The plan of this lecture is as follows. In Section [L2] I shall discuss possible acceleration mech- 



anisms of cosmic rays and astrophysical objects which potentially can be their sources. In Section 1.3 



I present the main energy loss processes for UHECR particles and briefly discuss their deflection in the 
magnetic fields. In Section [L4] I sum up recent observational results from the Pierre Auger Observatory 
and other experiments. In Section [T3] results on anisotropy at highest energy are discussed. In Sec- 
tion [L6] I review expectations on secondary photons and neutrinos from UHECR protons. Results are 



summed up in Section 1.7 



1.2 Acceleration 

There are several possible acceleration mechanisms that can work in astrophysical objects. These include 
first-order Fermi acceleration on the shocks in plasma or acceleration in the potential difference, which 
we call one-shot acceleration below. However, in any case, the Larmor radius of a particle does not 
exceed the accelerator size, otherwise the particle escapes from the accelerator and cannot gain energy 
further. This criterion is called the Hillas condition [6| and sets the limit 

£ <£n^ qBR (1) 

for the energy £ gained by a particle with charge q in the region of size R with the magnetic field B. 

The maximum energy of the accelerated particle can be restricted even more than required by 
Eq. ([T]) if one takes into account energy losses during acceleration. Unavoidable losses come from 
particle emission in the external magnetic field, which can be either synchrotron-dominated if the velocity 
of the particle is not parallel to the magnetic field, or curvature-dominated in the opposite case. 

In Fig. [2] in the plane magnetic field versus acceleration region size, the Hillas condition Eq. ([T]) 
is shown by a thick black line. The left figure is for protons and the right one for iron nuclei. Possible 
acceleration in different astrophysical objects is shown with thin solid figures. Notations are the follow- 
ing: NS are neutron stars, GRB are gamma-ray bursts, BH are black holes, AD are accretion disks, jets 
are jets in active galaxies, K and HS are knots and hot spots in the jets, L are lobes of radio galaxies, 
clusters are clusters of galaxies, starbursts are starburst galaxies, voids are voids in large-scale structure. 
Additional notations in brackets are subtypes of active galaxies: Sy for Seyfert galaxies, BL for BL Lac 
galaxies and RG for radio galaxies. Only objects above the Hillas line have the potential possibility to 
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Fig. 2: The Hillas plot with constraints from geometry and radiation losses for 10^^ eV protons (left) and iron 
(right). The thick line represents the lower boundary of the area allowed by the Hillas criterion, Eq. ([T]). Shaded 
areas are allowed by the radiation-loss constraints as well: light grey corresponds to one-shot acceleration in the 
curvature-dominated regime only; grey allows also for one- shot acceleration in the synchrotron-dominated regime; 
dark grey allows for both one-shot and diffusive (e.g., shock) acceleration. 



accelerate particles to 10^ eV. This is a necessary condition, but not enough for a specific acceleration 
mechanism. As seen from Fig. [2} for example, neutron stars cannot accelerate particles to highest en- 
ergies under any condition, while shock acceleration would work only for objects presented in the dark 
grey corner of this plot. 

1.3 Propagation 

Owing to expansion of the Universe, particles which come from sources at redshift z lose their energy as 

Ep^E'p^ Ep/{1 + z) . (2) 

A typical energy loss distance, i.e., distance at which particles lose a significant part of their energy for 
this process, is of the order of z ^ 1 (50% of energy), i.e., i? ^ 3 Gpc = 10^^ cm. 

As well as during propagation in the intergalactic space, protons lose energy due to two other main 
processes of interactions with Cosmic Microwave Background (CMB) photons. Those are electron- 
positron pair production and pion production. In both processes massive particles have to be produced 
and they have threshold energy. Since the typical energy of CMB photons is very small, ecMB = 
6 X 10~^ eV, the threshold for those processes is very high. Only at energies above Et^ = '^^/^cmb ^ 
10^^ eV does the electron-positron pair-production process become important: 

P + ^i^ + e+ + e-. (3) 

The typical energy loss distance for this process is 

Mp 1 97 

R = — ^ = 600 Mpc = 2 X 10^^ cm , (4) 

2me (Tpe+e-^CMB 

where ucmb = 400/cm^ is the density of CMB photons, crpe+e- ^ 10~^''/cm^ is the proton-pair 
production cross-section. The factor Mp/2me comes here from the fact that in every interaction a 
proton loses only a tiny fraction of its energy proportional to the proton/electron mass ratio. 
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Fig. 3: Left: Nucleon interaction length as function of energy from Ref. |7 |. Attenuation due to pion production 
Eq. ([5]) presented by the thick soHd Hne, the same for pair production Eq. ^ presented by the thin soHd Hne. 
Right: Horizon (maximal distance) from which protons with given or higher energy can arrive. Lines for 10%, 
30%, 50%, 70% and 90% of events are shown @. 



At energies above threshold Eth ^ ^ttMp/ ecMB = 10^^ eV, the pion production process domi- 
nates energy losses. This process for cosmic rays was first considered by Greizen, Zatsepin, and Kuzmin 
in 1966 |9| and is now named the GZK process. 

The typical energy loss distance for this process is 

R = ^ ^ = 100 Mpc = 3 X 10^^ cm , (6) 

where ap^ ^ 6 x 10~^^/cm^ is the proton pion production cross-section. The factor Mp/ttIt^ comes 
here from the fact that in every interaction the proton loses only 15-20% of its energy proportional to 
the proton/pion mass ratio. Note that at higher energies the dominating process in Eq. ([5]) is multi-pion 
production, in which the proton loses 50% of its energy in every interaction, however, the cross-section 
for this process aj-^^ = 10~^^/cm^ is a factor 6 lower than single pion production. 

None of the above processes allows a proton with high energy to come from a very large distance. 
The distance from which protons can come as a function of energy is presented in Fig. [3] (left) ||7|. 
The interaction length for pion production [Eq. ([5])] is shown by the dashed line. Attenuation due to 
pion production [Eq. ([6])] is presented by the thick solid line, the same for pair production [Eq. ([3])] is 
presented by the thin solid line. Figure |3] (left) shows the average distance travelled by a single particle. 
However, for searches of UHECR sources the important question is the maximum distance or horizon 
from which UHECR can come to the detector. In Fig. [3] (right) we present the horizon as a function of 
minimal proton energy. The lines 10%, 30%, 50%, 70% and 90% show the fraction of events which 
come from a given distance. For example, 90% of events with £^ > 10^^ eV should come from distances 
R < 100 Mpc. This distance is sometimes called the GZK distance, because energy losses in this case 
are dominated by the GZK process of Eq. ([5]). 

The dominant loss process for nuclei of energy E > 10^^ eV is photodisintegration A + j ^ 
(A — 1) + in the CMB and the infrared background due to the giant dipole resonance Ijlpl. The 
threshold for this reaction follows from the binding energy per nucleon, ~ lOMeV. Photo-disintegration 
leads to a suppression of the flux of nuclei above an energy that varies between 3 x 10^^ eV for He and 
8 X lO^*^ eV for Fe. 





Fig. 4: Sky map of the UHECR proton deflections for energy E = 40 EeV in two different models of Galactic 
magnetic field from Ref. |1H . The colours show deflections from to 10 degrees. 




Fig. 5: Fraction of the sky in which the deflection in the extra- Galactic magnetic field is bigger than the given 
value. Left: Constraint simulation of K. Dolag et al. p2| . Right: Simulation of Sigl et al. p3| . 



Since UHECR are charged particles, they not only lose energy in the interactions with background 
photons, but also when deflected by Galactic and intergalactic magnetic fields. 

The magnetic field of the Milky Way galaxy is conventionally modelled as a sum of the regular 
and turbulent components of the field in the disk and halo of the Galaxy. This means that the deflection 
in the Galactic field Ocai is a superposition of at least four terms: 



/J /jregular , /jturbulent , /jregular , /jturbulent 

C^Gal - C^Disk + ^Disk + ^Halo + ^Halo 



(7) 



The deflection angle of UHECR in a regular magnetic field after propagation of distance D is 
given by: 
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(8) 



where where B± is the magnetic field component orthogonal to the line of sight, -Euhecr is the particle 
energy, and Z is the atomic charge. In the case of deflection by the turbulent field on the distance D much 
larger than the correlation length of the field A b and where the deflection angle is small, the deflection is 
given by 
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(9) 



The deflection angles by regular and turbulent components of Galactic Disk and Halo, ^J^guiar 
and 6i°rt^uient ^^e given by Eqs. (jsj), Contributions of the Halo fields Of^^^^^^ and 9f^^{^^^^^^ are less 
certain, but the result in deflections is usually assumed to be less than the one for disk fields. 
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Deflections of UHECR by the regular field in the disk 6^^^^^^^ have been studied in many theo- 
retical models. A sky map of the deflections of UHECR with E = 40 EeV in two different models is 
presented in Fig. |4] Despite both models being consistent on average with the expectation of Eq. ([8]), 
predictions in any given direction are strongly model-dependent. 

Extragalactic magnetic fields are unknown except in the centres of galaxy clusters. Therefore one 
has to use theoretical models for the evolution of the magnetic fields. In such models magnetic fields 
follow the formation of large-scale structures. Because the structure of extragalactic magnetic fields is 
very non-trivial with very large fields near large-scale structures and tiny fields in the voids, one cannot 
use Eq. (|9]) everywhere. Instead one can introduce a fraction of the sky with deflections lower than a 
given value. Unfortunately, modern models give a very broad range of predictions. In Fig. [5] we show 
calculations by two groups that show very different results. The group of K. Dolag et al. made constraint 



simulations of local large-scale structures within 100 Mpc around the Earth |[T2|. This means that all 
big structures such as clusters of galaxies are located in exactly the same places as in the real sky. Also 
the density of points in this simulation is adaptive with more points at clusters and fewer on filaments. 
The results of this simulation are shown in Fig. |5] (left). According to this simulation, at 100 Mpc from 
the Earth only in 2% of the sky are deflections bigger than 0^^^^ = 1°. Those places are centres of 
galaxy clusters. In contrast the simulation of G. Sigl et al. | |T3| uses a uniform grid with more points 
on filaments and fewer on clusters. Unfortunately this simulation is not constrained and thus cannot be 
directly compared to local large-scale structures. In this simulation 0^^^^ > 50° in 60% of the sky. 

Let us note here that propagation energy losses are extremely important for the understanding of 
experimental results on spectrum and composition discussed in the next section. Deflections in turn are 



a key issue in the anisotropy studies in Section 1.5 



1.4 Observations 

In this section we shall discuss the present status of UHECR observations. 



Pilmsfv Ca^mlc Ray 



uv Fluoretcent pnotont 
Itotropic Emltfllor 




Charged partlctts 
£leclrcmagnetlc 
Shower 



Fig. 6: Extensive air shower produced by a UHECR particle in the atmosphere 



We start with the detection of UHECR in the atmosphere. The typical column density of the 
atmosphere is 1000 g/cm^ and in Ig there are N = lO'^^ protons, while a strong cross-section is app ^ 
10~^^ cm^. Thus UHECR protons or nuclei should interact within the atmosphere many times before 
they reach the Earth's surface. In these interactions it would produce extensive air showers. An example 
of such a shower is illustrated in Fig.|6] After first interaction the primary proton or nuclei would produce 
a large number of pions. Neutral pions would start an electromagnetic cascade, while charge pions would 
produce muons. At the maximum of shower development one expects N — 10^~^^ particles distributed 
in an area with a radius of a few kilometres. At this point the shower mostly consists of 10 MeV 
electrons and photons and only 5-10% of its energy is in muons. If the shower is not vertical, the 
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column density increases as 1/ cos{9 zenith) and reaches 2000 glcvc? for 9 zenith = 60°. At such a depth 
all electromagnetic components of the shower disappear and the shower would consist of muons only. 
Another way to detect the shower is to look for fluorescent UV light of nitrogen atoms in the atmosphere. 
This method is called 'calorimetric' and gives a 3-dimensional image of the shower. The main problem 
with this method is that detection is possible only on a moonless night, making the duty cycle of such 
detectors possible only 10-15% of the time. Finally one can detect direct Cherenkov light of the charged 
particles, but since this light is concentrated only within the central kilometre of the shower, one cannot 
use this technique at highest energies. 




Fig. 7: Pierre Auger Observatory detector with more than 1600 water tanks and 4 fluorescence telescopes (see 
Ref. (141 details) 



The Pierre Auger Observatory (Auger) is the largest UHECR detector in the world at the moment 
with an area of 3000 km^. Such a big area is required to collect enough statistics with UHECR at highest 
energies E > EeV, because the flux of such UHECR is tiny, one particle per 100 km^ per year. Auger 
is located on the plateau at an altitude of 1000 metres in the Mendoza province of Argentina. The ground 
detector consists of 1600 water tanks distributed 1.5 km from each other as presented in Fig. |7j Also 
there are four fluorescence telescopes pointed at the atmosphere above the ground detectors as shown 
in Fig. |7j Detection of 10% of showers both by fluorescence detectors (ED) and by ground detectors 
guarantees a good quality of events and at the same time allows one to calibrate the ground detectors by 
ED. 

In Fig. [8] we show a recent energy spectrum which was measured by Auger before 3 1 March 
2009 [4 |. The steeply falling flux of UHECR is multiplied by in order to show details of the spectrum. 
The total systematic energy error is 22% and is shown in the top right corner of the figure. For energy 
bins with £^ < 3 • 10^^ eV statistical errors are not important, while at highest energies E > EeV the 
shape of the spectrum is still uncertain and more statistics are needed. On the other hand, the suppression 
of the spectrum is statistically significant and is clearly seen in Fig. [8] 

This is an important experimental result, since it is independent confirmation of similar observa- 
tions made by the HiRes experiment [ 2 |. Thus cutoff in the energy spectrum exists. However, there are 
several questions to be answered before one can tell that this really is a GZK cutoff. First, is this cutoff 
due to the maximum energy of sources, or to energy losses? In Section 1.2 we have seen that indeed 
the maximum energy for many types of sources is close to 10^^ eV. The ultimate answer to this question 
would be the detection of several sources at different distances with cutoff following expectations of 
energy losses. 

Second, is the chemical composition of UHECR proton-dominated at those energies? Since in 
our Galaxy all elements up to iron are accelerated to energies around the knee E = 10^^ eV, the same 
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Fig. 8: Energy spectrum of UHECR as a function of energy measured by the Pierre Auger Observatory and model 
predictions for iron nuclei (blue) and protons (red) ||4|| 



situation can exist in astrophysical objects which accelerate to highest energies. Experimental answers to 
this question can be found in the future even by Auger, but already current data show that the composition 
becomes heavy at high energies. Indeed, recent Auger results from Refs. |[l5}[T6j are shown in Fig. [9) In 
Fig. [9] (left) we present the shape of the shower development in the atmosphere as seen by a fluorescence 
telescope. The signal is proportional to the number of electrons and positrons in the shower. Signals 
grow due to the development of electromagnetic cascades. The maximum of the signal corresponds to 
the maximum development of the cascade in the atmosphere. After that the shower loses its energy 
due to dissipation effects. The depth of the atmosphere corresponding to the maximum of the shower 
development is called X^^x- For the example presented, this maximum is at X^ax = 753 g/cm^ and the 
energy of the event is = 1.6 • 10^^ eV |j_5|. At the same energy, protons on average interact much 
deeper in the atmosphere than heavy nuclei. 

On the top panel of Fig. [9] (right) one can see the results of the most common hadronic models 
presented with red lines for protons and with blue lines for iron. The example event in Fig. [9] (left) is 
definitely proton-like. The averaged X^ax values in each bin are presented in the same figure for both 
the Auger and HiRes experiments. Both results are consistent with each other showing a relatively light 
composition from 10^^ eV to 10^^ eV. However, Auger shows heavier composition at highest energies. 

The main problem when measuring the composition with Xj^ax is its strong model dependence, 
as seen on the top panel of Fig. [9] (right). There are two complementary ways out. One is to use 
a composition-sensitive parameter that weakly depends on the model choice. Such a parameter is 



RMS{Xrnax) = V i^max) {^max)'^ ^ presented in the lower panel of Fig. [9| (right). One can see 
that according to this measurement the composition becomes heavier at high energy. Another impor- 
tant way is to test models and find the best one. For this purpose a dedicated experiment LHC forward 
(LHCf) was constructed at CERN. The idea of this experiment is to measure the neutral particles emitted 
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Fig. 9: Left: Measurement of shower development by signals in the fluorescence detectors as a function of depth 
in the atmosphere. The maximum of shower development in this example is X^ax = 753 g/cm^ p5| . Right: 
Average X^ax of showers measured by HiRes and Auger and RMS of X^ax measured by Auger in 2009 |T6| . 
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Fig. 10: The layout around the Interaction Point 1 (IPl) of the LHC. The structure at the centre indicates the 
ATLAS detector surrounding the collision point. The LHCf detectors are installed in the instrumentation slot of 
the TANs located ±140 m from IPl. Two independent detectors, LHCf Arml and LHCf Arm2 are installed at 
either side of IPl |[l7l. 



in the very forward region of LHC collisions at low luminosity. The configuration of this experiment is 
presented in Fig. [10] Data required for testing the hadronic models will be collected in the first scientific 
runs of the LHC 1 17 |. Thus in the near future we shall have better knowledge of hadronic models and 
more understanding of the composition of UHECR at highest energy. At present the Auger results indi- 
cate heavy composition; this was not confirmed by independent measurements and the fraction of light 
nuclei in the data remains uncertain. 

This is a very important question for searches of UHECR sources, which we shall discuss in the 
next section. 
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1.5 Anisotropy 

Since for every UHECR event the arrival direction is detected, for tens of years many attempts were made 
to find sources of UHECR in the experimental data. Unfortunately none of them has been confirmed so 
far. There are two ways to look for the sources. One is to look for the data itself and try to find anisotropy 
in autocorrelation factions. The second is to pick up a catalogue of possible sources and look for the 
cross-correlations with this catalogue. This second way always requires confirmation by an independent 
data set, since completeness of the catalogue is a very complicated issue and it is difficult to estimate the 
probability due to the parameter choice a posteriori. 



Here we start with autocorrelations. In the left panel of Fig. [TT] one can see the sky map with the 
arrival direction of events with E > 40 EeV in several old experiments, including SUGAR, AGASA, 
HiRes, Yakutsk, Havera Park, Volcano Ranch, and Fly's Eye. On the right panel of the same figure one 
can see the probability that autocorrelations between selected events are by chance within a given angle. 
One can see that the probability is minimal at angles 20-25 degrees. After penalization on the choice of 
angle, the probability that this happened by chance is P = 0.3% [ 18 1. This clustering of events on rather 
moderate scales can be due to the location of the sources in the Large Scale Structure. 



The same probability in the first Auger data is presented in Fig.[T2]as a function of both energy and 
angle. One can see that for exactly the same energy £^ = 40 EeV and angle 9 = 20°-25° the probability 
is P ^ 10~^. However, recent results with larger statistics did not show more significant anisotropy at 
such energies [20J . This makes the situation with anisotropy in the data less clear. 

Now let us discuss correlations with astrophysical objects. First Auger data have shown strong 
correlations with nearby active galaxies called Active Galactic Nuclei (AGN). Namely, 12 out of 14 
events with E > 57 EeV were correlated within < 3.1° from 472 AGNs from the Veron catalogue with 
distances R < 75 Mpc. This correlation was considered by the Auger Collaboration as a formal way to 
study the deviation of cosmic rays from isotropic distribution. Data from Period I was tested with the 
prescription during Period II, where 13 new events were detected, out of which 9 obeyed prescription 
parameters. The prescription was fulfilled, i.e., the observed sky was considered anisotropic at the 99% 
confidence level |21|. Data used in this publication and shown in Fig. [13] (left) correspond to Periods I 
(not shown) and II shown in Fig. [13] (right) before the vertical line. Unfortunately this correlation was 
not confirmed in the later data [Period III in Fig. [T3] (right)]. 

It does not mean that all anisotropy signals in Auger have completely disappeared. There is still a 
remaining excess of events around the Cen A galaxy on scales of 20 degrees, see Fig. [14] This anisotropy 
has to be tested by future data. 
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Fig. 12: Probability of autocorrelations as a function of energy and angular distance between events, see Ref. (19} 




Fig. 13: Left: Sky map of arrival directions of 27 UHECR with E > 57 EeV measured by the Pierre Auger 
Observatory before August 2007 in galactic coordinates (circles) and 472 nearby AGNs (red stars) |[2T|. Blue 
contours show the Auger exposure. Right: Likelihood ratio for events after formulation of the prescription. Period 
II is for data on the left panel. Period III is for new data up to March 2009 | [20| . 



1.6 Secondary photons and neutrinos from UHECR 

As was discussed in Section [13} protons lose their energy in pair production and pair production reac- 
tions. Since secondary pions quickly decay, secondary photons and neutrinos are produced. Neutrinos 
propagate to the Earth without interactions on the way, but photons cannot. They start to interact with 
background photons and produce pairs. Electrons and positrons in turn up-scatter CMB photons or pro- 
duce synchrotron radiation: 



T ~l~ ^background ^ ^ ~\~ G 

+ ^background ^6^+7 (10) 
+ B + synch 



The sequence of processes in Eq. ( [T0| ) is called an electromagnetic cascade. At energies above 
10^^ eV the cascade proceeds on the CMB background (400/cm^), but at lower energies pair production 
on CMB is impossible. At such energies the cascade continues on a much less abundant infrared back- 
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Fig. 15: Left: Fluxes of protons and secondary photons as a function of energy. Primary protons with spectrum 
XjE^-^ and maximum energy Emax = 10^^ are shown by the thin red Hne. Secondary protons fit the UHECR 
spectrum from E > 10^^ eV (thick red Hne). Secondary photons from all reactions are shown by the blue dashed 
line and from pion production only, Eq. Q by the magenta line. Right: Fluxes of UHECR and secondary photons 
in the case of iron nuclei primaries with spectrum 1/E'^-^ and maximum energy E^ax = 10^^ ^V. The remaining 
iron nuclei are shown by the green line. Secondary protons by the magenta line. Secondary photons by the blue 
line. 



ground (1/cm^) and at lower energies on optical background (0.01/cm^). Then it stops at the multi-GeV 
energies of gamma rays. 



In Fig. [T5] we plot primary cosmic-ray and secondary photon fluxes from primary protons (left) 
and iron (right) from Ref. [22 j . Secondary protons after interaction fit the UHECR spectrum from E > 



10^^ eV in Fig. |15| (left). Secondary photons cascade down to the GeV region. Only a small fraction of 
photons come from the pion production reaction (magenta dotted line). Most of the photons generated 
are from the e+e~ production reaction with total flux shown by the dash-dotted blue line. The number 
of secondary protons is much lower in the case of iron primaries, as shown by the dotted magenta line 
in Fig. [15] (right). As a result, the secondary photon flux in the GeV region is much smaller in this case, 
on the level of 0.2% of the EGRET measurement. Also very high energy photons are absent in this case 
due to low maximum proton energy. 

In Fig. [16] we compare the range of the electromagnetic cascade fluxes from UHECR with other 
possible astrophysical contributions in the EGRET band. Note that most of the uncertainty of the UHECR 
cascade flux comes from an unknown source evolution. The scatter for a given class of sources is thus 
much smaller, as seen from Fig. [T6]for the case of AGNs. 



In Fig. [17] (left) we plot the possible range of GZK gamma-ray fluxes for a given proton flux 
which fit the UHECR spectrum. The range of fluxes comes from the variation of possible values of the 
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Fig. 17: Left: Example of GZK photon flux from Ref. |23 1. UHECR protons fit the HiRes spectrum. Secondary 
neutrinos are shown by a green fine. The remaining secondary photons are in the range between the blue lines. 
Right: Experimental upper limits on the photon fraction in the UHECR spectrum from Ref. pO| . 



extragalactic magnetic field and the range of the models for extragalactic radio background. Also on the 
same figure the corresponding neutrino flux is shown by a green line. In Fig. [17] (right) we show the 
experimental upper limits on the fraction of photons in the UHECR flux. The range of possible GZK 
photon fluxes corresponds to protons with a range of power law injection spectra and source evolution 
fitting the UHECR spectrum. One can note that the current best upper limits of Auger are still above 
the range of expected theoretical values. On the other hand, existing limits already exclude some exotic 
models. 



1.7 Summary 

In the first lecture we briefly discussed many aspects of UHECR physics. 

Observed cosmic rays have energies up to 10^^ eV. Acceleration in astrophysical objects to such 
energies is a very non-trivial task and there are no objects in our Galaxy which can do this job. There are 
very few classes of exceptionally powerful objects in the Universe, some of which can be real sources of 
UHECR. Accelerated particles lose their energy in interactions with the CMB background and are also 
deflected by electromagnetic fields during their propagation from sources to the Earth. 
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There are three important experimental challenges in UHECR physics: the spectrum of cosmic 
rays, the chemical composition of cosmic rays, and the search for anisotropics in the sky with the ultimate 
goal of finding UHECR sources. 

The cutoff in the energy spectrum at highest energies E > 6- 10^^ eV has now been established 
by two independent experiments, HiRes and Auger. 

The most striking result of 2009 was evidence of heavy composition, shown by the Auger experi- 
ment at highest energies, Fig.|9] This result still needs independent confirmation. Also the interpretation 
of composition measurements is affected by uncertainty in the hadronic models. This question can be 
clarified in the near future by the LHCf experiment. 

Finally, most challenging is the search for UHECR sources. The last result in this direction was 
made by Auger in 2007. They found that the sky is anisotropic at the highest energies, at least at the 
99% C.L., by looking at the correlations with nearby AGNs. Unfortunately those correlations were not 
confirmed in the new data, and the only anisotropy excess remaining in the Auger data at the highest 
energies is an excess around the Cen A galaxy, see Fig. [T4j 

During energy losses the UHECR protons produce secondary photons and neutrinos. Most of the 
secondary photons cascade down to the GeV energies, where this contributes to the diffuse extragalactic 
background. An experimental search for the remaining gamma rays at highest energies E > 10^^ eV is 



challenging and existing upper limits are just above theoretical predictions, see Fig. [17 

Thus there are many unsolved problems in UHECR physics. They require both theoretical and 
experimental efforts in the near and more distant future. 

2 High-energy gamma rays 
2.1 Introduction 

In this lecture I shall discuss the theory of TeV gamma rays and recent observations made in this field. I 
shall give a brief introduction to the experimental detection techniques and present some selected results 
on the subject. For more detailed study I would like to recommend the recent review by F. Aharonian, 
J. Buckley, T. Kifune, and G. Sinnis 

Relativistic particles can travel with a speed larger than the speed of light in the medium V > 
Vm = c/n. Here n > 1 is the refractive index of the medium. This index in the air is na = 1.008 and in 
water nw 1.33. 

The charged particles polarize the molecules of the medium, which then return rapidly to their 
ground state, emitting prompt radiation called Cherenkov radiation. This radiation is emitted under a 
constant Cherenkov angle with the particle trajectory, given by 

X c 1 

cos5 = = = (11) 
V nV pn 

The minimal energy of a charged particle is 

_ Emin _ n 

Particles with higher energy will produce a cone of Cherenkov light. This effect is used by 
Cherenkov telescopes for air (H.E.S.S., MAGIC, Veritas, CTA) and by ground experiments in water 



(Milagro, HAWK). Detection of the shower in air and in water is illustrated in Fig. 18 



We present examples of such experiments in Fig. [19] On the left panel we show a view of the 
H.E.S.S. experiment. This experiment made the most significant contribution to the development of 
TeV gamma-ray astrophysics in recent years. On the right panel we show the Milagro experiment, a 
pioneering experiment in water Cherenkov techniques. 



14 





Fig. 18: Detection of high-energy gamma rays by Cherenkov telescopes in air (left) and in water (right) 




Fig. 19: Examples of gamma-ray experiments: Cherenkov telescope H.E.S.S. (left) and water pool Milagro (right) 
2.2 Point sources of TeV gamma rays 

TeV gamma-ray astrophysics is developing very quickly. One can see the number of detected sources in 



the sky as a function of time in Fig. [20] From 3 sources in 1995 one has 32 sources in 2005 and 80 sources 
in 2008. In addition not only does the number of observed sources grow, but also the number of different 
populations of sources. This is a very important fact for future experiments with better sensitivity like the 
Cherenkov Telescope Array (CTA). They would have a very large potential for detecting many different 
classes of sources. 



In particular, in Fig.|20]on the bottom panel, red circles show extragalactic sources which contain 
BL Lac objects, radio galaxies, and starburst galaxies. Also in the galactic plane there are many different 
classes of objects, which include supernova shells, pulsar wind nebulas, pulsars, binary systems and dark 
objects. Dark objects mean they were detected in gamma rays, but there is no corresponding source in 
other wavebands. 



The sensitivity of gamma-ray detectors to point sources as a function of energy is shown in Fig. 21 



The sensitivity of air telescopes is shown for 50 hours of observation for one source. The sensitivity for 
ground experiments is shown for 5 years, but they observe all the sky (27rsr). At low energies E < 
10 GeV the sensitivity of the GLAST (Fermi) satellite is the best. One year of observations are shown. 
At large energies £^ > 10 TeV the ground air- shower experiments (Tibet) have the best sensitivity. Future 
CTA projects will be orders of magnitude better than present-day experiments from 10 GeV to 10 TeV 
energies. 

Another important fact is that gamma rays cannot travel freely in the intergalactic space. They in- 
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1995 




Fig. 20: Sky in the TeV gamma rays with 3 sources in 1995 (top left), 32 sources in 2005 (top right), and 80 
sources in 2009 (bottom) 



teract with optical/infrared background photons and disappear producing pairs of electrons and positrons. 
In Fig. [22] one can see the main backgrounds for gamma-ray propagation. They are shown in units of 
photon density per cm^. The largest contribution comes from the CMB background with 400 photons 
per cm^. However, owing to the small energy of CMB photons, this background is important only for 
E > 1000 TeV. For the experimentally interesting energy range £^ < 100 TeV the main backgrounds are 
infrared and optical. Since those backgrounds are created by galaxies and partly by dust they are strongly 
model dependent both as a function of energy and as a function of redshift. 

Optical depth can be defined as 

r{E) = R ' ajj{E) • ntack{z, e) , (13) 

where R is the distance travelled by photons, a^^{E) is the pair-production cross section, and ni^acki^i ^) 
is the density of background photons. Distances on the cosmological scale are often expressed in terms 
of redshift. One can express it through the Hubble law R = z • c/Hq, where Hq = 70 km/s/Mpc is 



the Hubble constant. In Fig. 23 contours of constant optical depth r{E) are shown on the plane redshift 
versus energy for r{E) = 1, 3, 10 in two different models of IR/O background. 

There is one important difference between air Cherenkov telescopes and water Cherenkov detec- 



tors. In Fig. 24 we plot world-wide monitoring of the nearby BL Lac object Mkn 421 as a function of 
time. One can see that air Cherenkov telescopes can see a signal only on moonless nights, which restricts 
their operation to the corresponding intervals of time. On the contrary, water Cherenkov telescopes op- 
erate all the time they can see a source, which will allow source activity to be detected all the time. 
On the other hand, the problem of water Cherenkov experiments is poor sensitivity, which will prevent 
them from detection of relatively low fluxes and very fast variations in time. Thus both techniques are 
complementary to each other. 
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Fig. 21: Sensitivity of gamma-ray detectors to point sources, from Ref. (24] 
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Fig. 22: Redshift for gamma rays as a function of energy. Lines show constant optical depth in two models of 
IR/0 background. 



In Fig. |25]one can see a view of the central part of the Milky Way galaxy in three energy bands: 
optical, infrared, and TeV gamma rays. At least three astronomical source populations: supernova rem- 
nants (SNRs), pulsar wind nebulae (PWNe), and binary systems (BSs) are represented in this figure. 
In addition, the H.E.S.S. observations of the central region of our Galaxy revealed a diffuse TeV 7- 
ray emission component which is apparently dominated by contributions from giant molecular clouds 
(GMCs). These massive complexes of gas and dust most likely serve as effective targets for interactions 
of relativistic particles from nearby active or recent accelerators. Thus one may claim that four galac- 
tic source populations are already firmly established as effective TeV 7-ray emitters. Meanwhile, many 
sources discovered by H.E.S.S. in the galactic plane remain unidentified. Although some of these sources 
might have direct or indirect links to SNRs, PWNe, and GMCs, one cannot exclude that a fraction of the 
H.E.S.S. unidentified sources are related to other source classes. 

The Milagro telescope has made the first measurement of the diffuse TeV gamma-ray flux from the 
Galactic Disk. Figure |26] shows the Galaxy (as visible from the Northern Hemisphere) in TeV gamma 
rays. In addition to the individual sources discussed above, the image (compiled from Milagro data) 
shows the existence of a diffuse TeV gamma-ray flux between galactic longitudes of 30° and 90°. 
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Fig. 23: Redshift for gamma rays as a function of energy. Lines show constant optical depth in two models of 
IR/0 background. 
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Fig. 24: Observation of Mkn 421 as a function of time 



2.3 Extragalactic magnetic fields 

Another very important field which will benefit in the near future from TeV gamma rays is Extragalactic 
Magnetic Fields. 

Indeed, as discussed above, TeV gamma rays emitted by astrophysical sources can be measured by 
detectors on Earth. Practically all TeV gamma rays from galactic sources come directly to the detectors. 
However, this is not true for extragalactic sources. As one can see from Fig.[23j even for nearby sources 
like Mkn 501, gamma rays with £^ > 10 TeV cannot come freely to the detector. The pair production on 
Extragalactic Background Light (EBL) reduces the flux of 7-rays from the source by 



(14) 



where F{Ejq) is the detected spectrum, Fq{E'^^) is the initial spectrum of the source, and t{Ejq^ze) 
the optical depth Eq. ( [T3] ). The typical distance which a primary gamma ray travels is 

-1 



is 



(1 + ^) 



20 TeV 



Mpc , 



(15) 



where a numerical factor k = k{E. 



70 ' - 



1 accounts for the model uncertainties. 
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Fig. 25: Central part of the Milky Way galaxy in infrared, optical, and in TeV gamma rays. The TeV gamma-ray 
sky from H.E.S.S. observations with a large number of sources. 
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Fig. 26: The Milky Way galaxy in TeV gamma rays from galactic longitude 20° to 220° and galactic latitude from 
— 10° to 10°. The image is the culmination of a seven-year exposure by the Milagro instrument. 



The cascade electrons lose their energy via Inverse Compton (IC) scattering of the CMB photons 
within the distance 



2^3 



10^^(1 + 



77/ 



10 TeV 



cm 



(16) 



The deflection angle of the e+e~ pairs, accumulated over the cooling distance, depends on the correlation 
length of the magnetic field, A^- Note also that electrons and positrons travel much shorter distances than 
primary photons. 
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D 

Fig. 27: Detection of EGMF through observation of secondary emissions around a point source p5| 



The e+e~ pairs produced in interactions of multi-TeV 7-rays with EBL photons produce sec- 
ondary 7-rays via IC scattering of the Cosmic Microwave Background (CMB) photons. Typical energies 
of the IC photons reaching the Earth are 



11) 



771/2 



0.32 



E' 



70 



20TeV 



TeV 



(17) 



where ^cmb = 6 x 10~'^(1 + z^^) eV is the typical energy of CMB photons. In the above equation we 
have assumed that the energy of a primary 7-ray is E'^^ 2^ 2E'^ with E'^^ being the energy of the primary 
7-rays at the redshift of the pair production. Upscattering of the infrared/optical background photons 
gives a sub-dominant contribution to the IC scattering spectrum because the energy density of CMB is 
much higher than the density of the infrared/optical background. 

Deflections of e+e~ pairs produced by the 7-rays, which were initially emitted slightly away from 
the observer, could lead to 'redirection' of the secondary cascade photons toward the observer. This effect 
leads to the appearance of two potentially observable effects: extended emission around an initially point 
source of 7-rays |[25l-[27l and delayed 'echo' of 7-ray flares of extragalactic sources P28ll29fl. 




Fig. 28: Model predictions and estimates for the EGMF strength. Cyan shaded region excluded by present day 
measurements. Black ellipses show measurements of the field in the Galaxy and galaxy clusters. Left panel: left 
and right hatched regions show theoretically allowed range of values of (Ab,B) for non-helical and helical fields 
generated at the epoch of electroweak phase transition during radiation-dominated era. Middle panel: left and 
right hatched region show ranges of possible (Ab,B) for non-helical and helical magnetic fields produced during 
the QCD phase transition. Right panel: hatched region is the range of possible (Ab,B) for EGMF generated 
during recombination epoch. Dark grey shaded region shows the range of (Ab,B) parameter space accessible for 
the 7-ray measurements via 7-ray observations pQ| . 

The above processes are illustrated in Fig. [27j Electron deflection 5 depends on the magnetic field 
in the region of deflection. Note, that, in principle, EGMF depends on the redshift, B' — B'{z). In the 
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simplest case, when the magnetic field strength changes only as a result of expansion of the Universe, 
B'{z) ~ i?o(l + z)"^, where Bq is the present epoch EGMF strength. This gives 
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Knowing the deflection angle of electrons, one can readily find the angular extension of the sec- 
ondary IC emission from the e+e" pairs 
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This is a key point for detection of the field, since extended emission depends on energy in a well-defined 
way and can be reconstructed using independent measurements at different energies. 



The possible ranges of the (A^,B) parameter space are shown in Fig. |28]for the cases when mag- 
netogenesis proceeds during electroweak or QCD phase transitions or at the moment of recombination. 

It is interesting to note that predictions for the strength and correlation length of the primordial 
magnetic fields fall in a region of (A^,B) parameter space which is not accessible for the existing mea- 
surement techniques, such as Faraday rotation or Zeeman splitting methods. However, it turns out that 
this region of (A^,B) parameter space is accessible for the measurement techniques which exploit the 
potential of the newly opened field of very-high-energy (VHE) 7-ray astronomy pO|. 



2.4 Summary 

Gamma-ray astronomy works, hundreds of sources have been detected in the GeV energy range and 
about one hundred in TeV energies. 

There are several major questions to be answered in the near future: 



- One needs to understand the hadronic component in a variety of astrophysical sources. 

- Extragalactic IR/O backgrounds have already been constrained by observations of TeV sources 
to factor two uncertainty. The next step is precision determination of those backgrounds using 
measurements of many sources at different redshifts. 

- For the first time one has a possibility to study primordial magnetic fields through TeV gamma-ray 
measurements. We can test models of primordial magnetic fields in the near future. 

There are several other important issues which were not discussed in this Lecture due to lack of 
time. The corresponding questions are: 

- Good measurements of blazar flairs can help to understand gravity near black holes. 

- TeV gamma rays give one more constraint/signature on Dark Matter. 

- Constraints on exotic physics (LIV, etc.) will be improved. 
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3 High-energy neutrinos 
3.1 Introduction 

In this lecture we discuss theoretical predictions and experimental efforts to detect Ultra-High Energy 
neutrinos. In Section |3.2| we discuss possible ways to detect UHE neutrinos and their corresponding 
experiments. In Section [33] we show theoretical predictions for UHE neutrino fluxes and present the 
status of experimental searches for such fluxes. In Section we discuss another possibility to detect 
Galactic neutrino sources at multi-TeV energies. In Section |3.5| we summarize all the results of this 
lecture. 



3.2 High-energy neutrino experiments 

There are three types of ultra-high energy (UHE) neutrino experiments. 

First, neutrinos can be detected by UHECR experiments. There are two possibilities for this. First, 
one can use the fact that the atmosphere horizontally has depth 36 times the vertical depth. Relatively 
young electromagnetic horizontal showers can be caused by neutrinos only. Hadronic showers at such a 
depth consist of muons only. Second, one can look for events penetrating the Earth in the tau-neutrino 
channel, i.e., look for upward-going events. This was used by the Auger experiment (see Fig. [7]). The 
resulting limit on neutrino flux is shown in Fig. [31] Also less significant limits were presented by previous 
UHECR experiments including Fly's Eye, AGASA, and HiRes (the HiRes limit is also shown in Fig. [31]). 




Fig. 29: IceCube detector. Left: Configuration of the IceCube detector. Eighty strings will be located at a depth 
of 1.5 km in the Antarctic ice filling a volume of one cubic kilometre. The present construction stage is also 
shown pT| . Right: Simulation of a high-energy neutrino event in the IceCube detector |32|. 



Second, one can detect neutrinos in the water or in the ice by detecting Cherenkov light created by 
corresponding leptons after neutrino interaction in the medium. There are two important backgrounds 
for such measurement. First, secondary leptons, mostly muons, should not be confused with secondary 
muons from extensive air showers in the atmosphere. In order to reduce the background of atmospheric 
muons one has to put the detector at a depth greater than one kilometre from the surface. Second, there 
are atmospheric neutrinos created by the same cosmic rays, which would produce isotropy in the space 
energy-dependent background. In order to fight this background, one either has to go to high energies 
E > 10^^~^^ eV, where it is small, or look for point sources on top of this background. 

Experiments that worked with these techniques in the past were Baikal and ANTARES in water 
and AMANDA in ice. All those experiments had a volume 0.1 km^ or less. The new-generation ex- 



22 



periment IceCube with a volume of 1 km^ is in the construction stage at the moment. In Fig. 29 in the 
left panel one can see the configuration of this experiment, which consists of 80 strings, filling a cubic 
kilometre volume in the Antarctic ice at a depth of 1.5 km from the surface. Strings already implemented 
are shaded blue on top of the picture (see Ref. pT[ for more details). Also, as shown in the figure the 
top of the detector is covered by an array of ice tanks (ice top). In the right panel one can see a Monte 
Carlo simulation of a high-energy neutrino event, detected by the IceCube experiment. First results of 



this experiment will be discussed in Section 3.4 




/ 

Cover (partially cut away) 



Fig. 30: ANtarctic Impulsive Transient Array (ANITA) radio balloon experiment. Array of radio antennas flying 
in the ballon, as shown on the left panel. It flies in circles over the Antarctic ice at a height of 37 km (see right 
panel) and looks for radio signals which UHE neutrinos create in the ice. 



Finally, radio neutrino experiments exploit the Askaryan effect in which strong coherent radio 
emission arises from electromagnetic showers in any dielectric medium. High-energy neutrinos trigger 
a cascade of electromagnetic particles in the medium, which has net charge and can emit an analogue 
of Cherenkov light in the radio energy range. The main point of this effect is that the length of the 
radio wave is macroscopic (tens of centimetres) and is bigger than the size of the cascade itself. This 
in turn means that all electrons in the cascade emit coherently. The effect was first observed in 2000 at 
SLAC. Recently the Askaryan effect has been clearly confirmed and characterized for ice as the medium, 
as part of the pre-flight calibration of the ANITA- 1 payload. The Askaryan effect can be seen only at 
high energies E > 10^^~^^ eV. Experiments using this effect benefit from the absence of atmospheric 
neutrino flux at such high energies, but they also have to look over a huge effective volume in order to 
see tiny neutrino fluxes at highest energies. 



Experiments that used this effect to search for UHE neutrinos are FORTE | [33| , RICE | |34| , and 
ANITA | [35l[36| . FORTE is a satellite experiment, which, in particular, looked over the Greenland ice. 
Unfortunately, the threshold of this experiment was very high, Ej^ > 10^^ eV, so it could test only exotic 
top-down models. RICE was an array of radio antennas located in the ice at the South Pole at the same 
place as the AMANDA experiment. This experiment presented its final results in 2006. Finally, the 
most advanced for the moment of this kind of experiment is the ANtarctic Impulsive Transient Array 



(ANITA) radio balloon experiment, see Fig. |30j In the left panel one can see an array of radio antennas 
in the balloon. In the right panel one can see a schematic map of flight over the Antarctic at a height of 
37 km. 



3.3 Search for cosmogenic neutrinos 



As discussed in Section 1.3 UHECR protons lose their energy in interactions with CMB photons and 
produce pions at energies above threshold E > 6 • 10^^ eV. This GZK threshold was found in 1966 ||9). 
As long ago as 1969 Berezinsky and Zatsepin suggested that one can try to observe secondary neutrinos 
from pion decays and called them cosmogenic neutrinos p7| . Recently the ANITA Collaboration pro- 
posed to call such neutrinos Berezinsky-Zatsepin neutrinos, or BZ neutrinos [ |36| |. Below we follow this 
suggestion. 
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Fig. 31: Predictions of cosmogenic neutrino fluxes and theoretical bounds on them |[38l[39| 



One can calculate the flux of BZ neutrinos theoretically, after fitting the corresponding proton 
spectrum to the experimental flux above some energy. The absolute limit for neutrino flux comes from the 
fact that gamma rays unavoidably produced from tt^ decays and from electrons tt^ decays cascade down 
to GeV energies and the maximum flux of such gamma rays cannot overshoot the EGRET measuremen 



shown in Fig. 16 This bound on the BZ neutrino flux is called "gamma-ray bound" in Fig. 31 Note 
that there are many additional ways to create photons in the EGRET energy range, including electron- 
positron pair production discussed in the previous section, so the real BZ neutrino flux is always lower 
than this region. 
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Fig. 32: Experimental limits on cosmogenic neutrino flux. Best up-to-date ANITA- 1 limits based on no surviv- 
ing candidates for 18 days of live time shown as ANITA-2008 |36|. Also limits from Auger |4Q| , HiRes j41|, 
FORTE |^33J, Anita prototype ANITAlite |35 1, RICE and AMANDA II |,42J are shown. 



Also in Fig. [31] we plot two theoretical limits derived under a set of theoretical assumptions. One 
is called the Waxman-Bahcall (WB) bound and the other the MPR bound. On the same figure we show 
several examples of theoretical neutrino fluxes which violate both WB and MPR bounds, but all of them 
are consistent with the experimental gamma-ray bound. 
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In Fig. 32 we show present-day experimental bounds confronting theoretical predictions for BZ 
neutrinos from Ref. |36|. One can see that the best up-to-date experimental bounds come from the 
ANITA experiment. ANITA-1 was able to view a volume of ice of ^ 1.6 Mkm^ during 17.3 days, how- 
ever, volumetric acceptance to a diffuse neutrino flux, accounting for the small solid angle of acceptance 
for any given volume element, is several hundred km^ water-equivalent steradians at Ejy = 10^^ eV. This 
allowed them for the first time a tough theoretically interesting region, excluding part of the parameter 
space with highest neutrino fluxes. 

On the same figure one can see existing limits on diffuse neutrino flux from the Auger | [40| , 
HiRes ||4T|, FORTE |33|, Anita prototype ANITAHte |[35|, RICE |[34j, and AMANDA II (4T\ experi- 
ments. 



Let us note also that in Fig. [32] the composition is assumed to be proton-dominated. If recent 
Auger results presented in Fig. [9] are confirmed, theoretical expectations for neutrino flux in Fig. 32 will 
be strongly reduced. This will make observations of the diffused flux of UHE neutrinos an even more 
complicated issue. However, at lower energies one still can have a hope of seeing point sources with 
neutrinos, as will be discussed in the next section. 



3.4 Point sources of UHE neutrinos 
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Fig. 33: Neutrino-nucleon cross- section as a function of the neutrino energy. Charge-current and neutral-current 
contributions to the cross-section are shown with thin solid and dashed lines. The total cross-section is presented 
by a thick solid line. See Ref. [43j for details. 



At highest energies the neutrino flux is too low to detect one single source of neutrinos, but at lower 
energies E < 1000 TeV the flux from a single source can be high enough to detect it. Indeed, in Fig. [33 
the neutrino-nucleon cross-section is shown as a function of energy. This cross-section is proportional to 
E at low energies E < 1 TeV and to E^-^ at high energies E > 10^ GeV. Good candidates for neutrino 
sources in the Galaxy are objects emitting TeV gamma rays. They can produce neutrinos in the proton- 
proton collisions in objects in the case of binary systems and in the interaction with molecular clouds in 
the Galaxy. In the 10 TeV energy range 



ap^(lOTeV) = 10" 



cm^ 



(20) 



In the IceCube detector only a small fraction of neutrinos will produce a signal: 



(21) 
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where njcE ^ 10^^ /cm^ is the density of the ice and i? = 1 km is the height of the IceCube detector. 
The expected flux of neutrinos produced in the proton-proton collisions in the Galactic sources is 

F^^F^ = 10-^2^^^ ^ 3 . 10^^^ . (22) 
cm^s km yr 

Thus in the IceCube detector one can expect three events per year for a 10 TeV neutrino flux. 




Fig. 34: Left: Simulated detection of Milagro TeV galactic sources by IceCube. Right: Significance of Milagro 
hotspots after five years of observation of IceCube. 




Fig. 35: Equatorial sky-map of events (points) and pre-trial significances (/7-value) of the all-sky point source 
search in the 22-string IceCube detector |44|. The solid curve is the galactic plane. The most significant spot 
arrives in a random sky with probability P ~ 1%. 



We now present recent results for point-source searches using data recorded during 2007-08 with 
22 strings of IceCube (1/4 of the detector). An all- sky search within the declination range —5° to +85° 
found the most significant deviation from the background at 153.4° r.a., 11.4° dec. Accounting for 
all trials in the point-source search, the final value for this result is 1.34%, consistent with the null 
hypothesis of background-only events at the 2.2cr level. No obvious source candidates are near this 
location, and an analysis of the timing of the events did not find any evidence of a burst in time. The 
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location can be added to the a priori source candidate list for analysis using future IceCube data, in which 
case a similar excess would be identified with much higher significance [ |44| . 

3.5 Summary 

IceCube is half-complete. If it observes first sources, a new field of astroparticle physics will be started: 
neutrino astrophysics. If not, much bigger detectors are needed with a size of at least 10 km^. Secondary 
neutrino flux from UHECR protons can be detected by future radio experiments, like ANITA. Neutrinos 
from some bright galactic sources can be detected by IceCube. Extragalactic sources can be observed 
during bright flair activity. In order to detect continuous flux from sources like Cen A one needs detectors 
much larger than 1 km^. Galactic SN can be detected with neutrinos at low and high energies. Cubic- 
kilometre water detectors will be constructed if IceCube gives positive results. 
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